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Effects of Quenching and Partial Quenching on QCD Penguin Matrix 
Elements 
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We point out that chiral transformation properties of penguin operators change in the transition from un- 
quenched to (partially) quenched QCD. The way in which this affects the lattice determination of weak matrix 
elements can be understood in the framework of (partially) quenched chiral perturbation theory. 



1. (Partially) Quenched QCD penguins 

We consider AS = 1 LR operators of the form 



(1) 



with q = u,d,s, {qiq2)L,R = qil^PLM2, and 
Pl,r = (1 T 75)/2 left- and right-handed projec- 
tors. Color indices can be contracted in two ways, 
corresponding to the QCD penguins Q^fi. 

The operators in eq. (|l|) are obtained by the 
unquenched QCD evolution of the weak oper- 
ator from the weak scale ^ Mw down to the 
hadronic scale ~ rric- At hadronic scales, ma- 
trix elements of these four-quark operators are 
usually computed on the lattice in quenched or 
partially quenched QCD ((P)QQCD). The key 
point we wish to address here is that the chiral 
transformation properties of QCD penguin oper- 
ators change in the transition from unquenched 
to (P)QQCD. This has non-trivial consequences 
Q which can be understood in the framework 
of (partially) quenched chiral perturbation the- 
ory ((P)QChPT) §. 

PQQCD can be systematically formulated in a 
lagrangian framework by coupling the gluons to 
three sets of quarks ||^ : K valence quarks qm with 
masses m^i, my2i ■ ■ ■ , frivK, N sea quarks qsi with 
masses nisi, ^1^2, ■ ■ ■ , itIsn, and K ghost quarks 
qgi with masses mvi,my2, ■ ■ ■ ,nijjK- The ghost 
quarks are identical to the valence quarks, except 
for their statistics, which is chosen to be bosonic 



10]. Quenched QCD corresponds to the special 
case N = (no sea quarks). It can be shown 
that unquenched QCD below the charm threshold 
corresponds to the choice K = N = 3 and nisi = 
m^,, i = 1, . . . ,3 §,|]. 

The total number of quarks in PQQCD is thus 
2K + N, and correspondingly, the chiral symme- 
try group enlarges from SU{3)l x SU{3)ii to the 
graded group SU{K + N\K)l x SU{K + N\K)r. 
As a consequence, the operators of eq. (|l|), which 
are singlets under SU{3)r, are no longer singlets 
under the enlarged group SU{K + N\K)b.- They 
can be decomposed as 
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, / K K K 
A = diagll ,...,1 , ,. 

where the first K (valence) entries of A are equal 
to 1 — K/N, and the next N + K (sea and 
ghost) entries are equal to —K/N . The super- 
scripts PQS and PQA indicate that these oper- 
ators transform in the singlet and adjoint repre- 
sentations of SU{K + N\K)ii, respectively. In 
the quenched case the situation is special. The 
decomposition reads 



^QCD 
^penguin 
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N = 



2 ^penguin 

idiag(l, . . . 



^penguin ' 



1,-1, 



(4) 

,-1), (5) 



where the first K (valence) entries of N are equal 
to ^, and the last K (ghost) entries are equal to 
— ^. The first operator in the decomposition is a 
singlet, while the second is not, under SU{K\K)ii 
(NS for non-singlet). However, the unit ma- 
trix has now a vanishing supertrace, while N has 
not, and Q^^^,^ can mix with Q%g^,^ through 
penguin-like diagrams. 

2. Representation of QCD penguins in 
PQChPT 

QCD penguins transform as (8, 1) under the 
chiral symmetry group and start at order in 
ordinary ChPT (see e.g. ref. Q). Denoting the 
adjoint representation of the PQ group by A, 

Pw S 

we found in the previous section that Qpenguin 
transforms as (A, 1), while Qp^iguin transforms as 
(A, A) under SU{K + N\K)l ^ SU{K + N\K)r. 
To lowest order in PQChPT and in euclidean 
space, these operators are represented by 



,PQS 
penguin 

,PQA 
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(S 1) 

■ ~al ' str (AL^Lfj^) + a2 



(6) 



where = iT.dpT.\ X± = 2Bo{T,M'' ± MEt)^ 
with M the quark-mass matrix, Bq the param- 
eter Bq of ref. 0, E = exp(2i$//) describing 
the partially-quenched Goldstone-meson multi- 
plet, and / the bare pion-decay constant normal- 
ized such that /jr = 132 MeV. The low-energy 
constants (LECs) a^i2^ also appear in the un- 
quenched theory. However, the LEG a^^'*^ only 
appears in the PQ case, and multiplies an order 
p° operator, which is possible because A is non- 
trivial. It has a direct physical meaning, because 
the same LEC also appears in the bosonization of 
the EM penguin (which belongs to the same ir- 
reducible representation as Qp2i^uin)- Quenched 
bosonization rules to leading order are given by 
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'i^^^str {KL^L^) + a^f ^str {kX+) 
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The quenched case differs from the partially 
quenched one in two ways. First, all ChPT LECs 
depend on N, and thus their quenched values are 
not necessarily equal to their N values. Sec- 
ond, in this case the operator Qp^g^i^ and the 
EM penguin do not belong to the same irrep and 
their LECs are in principle not related jl]. 

3. Kaon matrix elements in PQChPT 

Since the new operators Q^en^um in eq. (||) and 
QpeZgmn in eq. (0) are of order / in ChPT (while 
singlet LR operators are of order p^), they po- 
tentially lead to an enhancement relative to the 
unquenched case. As the simplest example, we 
consider if — > tt and i^T — > matrix elements, to 
leading order in ChPT. New contributions also 
show up in direct K tttt matrix elements |^. 

It turns out that Q^^,^^,^ (and QpZ^mn) do 
not contribute to any of these matrix elements 
(i.e. K TT^ K ^ Q and K tttt) at or- 
der but they do in general contribute at or- 
der Since Qp^g^in starts at order p^, the 
new contributions from Qp^^iguin compete at the 
leading order of the chiral expansion of these ma- 
trix elements, and have to be taken into account 
even if one analyzes lattice results using only 
leading-order ChPT. For the PQ X — > tt matrix 
element, with degenerate valence quark masses 
(A//2 = M]^ = Af2 ^ 2Bomy), we find at order 



ipenguin 
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s,i) 



.,(8^1) 



where the LECs /3i 2 3 appear in the bosonization 

of Qp2i^uin at order [|l| . The A' ^ matrix 
element with arbitrary valence and sea quarks is 
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Here L{M) — log-^, and the result is given in 
the M S seheme, with A the running scale. M^si 
{M2si) is the mass of a meson made out of the 3rd 
(2nd) valence {i.e. the strange (down)) quark and 
the ith sea quark; analogously for {M2vi) 
with sea replaced by valence. At the order we are 
working, 



3si 



ML 



ML 



M^, which follows from M^^^ 



ML 



Ml 



Boiniys + iTisi), 
etc. These results also contain the unquenched 
result as a particular case, i.e. for iV = if = 3 
and by equating sea and valence quark masses, 
nisi = rrim, i — 1,...,K. Quenched {N — 0) 



results are obtained by replacing a\ 



(8,1) 



and /3{' 



(8,8) 



91,2 
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and by dropping all terms containing sea quarks. 

We conclude that (in general) the bosonization 
of the new operator Qp2^uin contributes to a 
weak matrix element at leading order in ChPT 
with non-analytic terms of the form KP log Af ^ , 
which are absent in the unquenched case, and 
with analytic contributions parameterized by new 
LECs /35*i [N ^ 0) or /3,^|,3 (iV = 0). 



4. Lattice strategies 

To leading order in ChPT, the physical K — > 
TTTT amplitude is determined by af'^\ and there- 
fore this is the LEG one wishes to extract from 
the K —> TT matrix element. It is then clear that 
this can be done by considering only the singlet 
penguin, Q^Zlmn^ in the PQ theory with TV = 3 
light sea quarks. This is equivalent to omitting 
all Wick contractions in which q and q in eq. (|^) 
are contracted (eye graphs), except when they 
correspond to sea quarks.0 Of course, since the 
K TT matrix element only gives the linear com- 
bination q;^*'^^ — aj^'^"* , the if — > matrix element 

of Qpe«3«m is also needed, as usual g. 

In the case that TV 7^ 3, one can still follow 



^The remaining "eye graphs" are those for which q (q) is 
contracted with s (d) in eq. (h 



a similar strategy in order to determine af'^\ 
However, since the LECs depend on the number 
of light sea quarks N, there is no reason why the 
result should be the same as that of the real world, 
which has = 3. In particular, in the quenched 
case different strategies are possible. If one as- 



sumes that a[^'^\N = 0) « a\^'''{N = 3), one 
should use the strategy described above. Alter- 
natively, one may include both the singlet and 
non-singlet operators, also in the conversion to 
K TTTT Ql . Since it is not known which (if any) 
of these strategies is a better approximation, we 
believe that the difference between them should 
be interpreted as an indication of the quenching 
error. 

Finally, we stress that our observations also ap- 
plies to LL penguin operators (in which case the 
new effects start at order fl],!!) and, more in 
general, to any weak matrix element to which 
penguin operators contribute, such as also non- 
leptonic B decays. 

We would like to thank the RBC collaboration, 
Claude Bernard, Guido Martinelli, Santi Peris 
and Steve Sharpe for discussions. 

REFERENCES 

1. M. Golterman and E. Pallante, 



(8.1) - 



hep-lat/0108010 . 

M. Golterman and E. Pallante, in progress. 
G. Bernard and M. Golterman, Phys. Rev. 
D49 (1994) 486. 

A. Morel, J. Physique 48 (1987) 111. 

S. Sharpe and N. Shoresh, Phys. Rev. D62 



(2000) 094503; |hep-lat/0108003 . 
J. Donoghue, E. Golowich and B. Holstein, 
"Dynamics of the Standard Model" (Gam- 
bridge, 1992); G. Bernard, in "From Actions 
to Answers, " proceedings of TASI '89, eds. 
T. DeGrand and D. Toussaint (World Scien- 
tific, 1990). 

J. Gasser and H. Leutwyler, Nucl. Phys. B25G 
(1985) 465. 

G. Bernard, T. Draper, A. Soni, H. Politzer 
and M. Wise, Phys. Rev. D32 (1985) 2343. 



